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Abstract. The tomographic method is employed to investigate the presence of 
quantum correlations in two classes of parameter-dependent states of two qutrits. 
The violation of some Bell's inequalities in a wide domain of the parameter space 
is shown. A comparison between the tomographic approach and a recent method 
elaborated by Wu, Poulsen, and M0lmer shows the better adequacy of the former 
method with respect to the latter one. 



1. Introduction 

In quantum mechanics, the correlation of physical observables can obey some 
constraints expressed in terms of inequalities. For example, the correlations of 
position and momentum of a particle provide the lower bound for the product of 
the uncertainties of these two quantities as given by the Schrodinger-Robertson 
uncertainty relation [1] [5] . 

For discrete variables like spins, the quantum correlations for entangled states, 
e.g., of two qubits, also provide some inequalities known as Cirel'son 's bound 2^2, 
which differs for similar CHSH-Bell bound 2 J^jSj associated with classical correlations 
typical for separable states of two qubits. 

In principle, the existence of bounds of this kind can be checked experimentally. 
The violation of Bell's inequalities, that is the violation of the bound of Bell's number 

2, was indeed proved experimentally (see, for example, ^j). 

Statistical characteristics as means, dispersions, correlations for quantum 
observables can easily be formulated in the probability representation of Quantum 
Mechanics (see, recent review [7]). This representation for spins was introduced 
in [SI inj- The quantum states in the probability representation are described by 
fair measurable probability distributions (called tomograms), which contain complete 
information on density operators of the states. In view of the possibility to use the 
standard formulae of classical probability theory in the probability representation for 
quantum-mechanical calculations with tomograms, there might exist some advantages 
to better understand the role of quantum correlations and of the corresponding 
quantum bounds for separable and entangled states, if one uses the tomographic- 
probability approach. 

Our goal in this paper is to investigate the existence of nonclassical correlations in 
selected two-qutrits states through Bell's number B using the tomographic approach. 
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2. A short review of the quantum tomographic method 

Quantum states are usually described in the canonical quantum picture by vectors 
in a complex Hilbert space or by density operators, describing pure and mixed states 
respectively. Together with the usual Schrodinger and Heisenberg picture there exist 
another possible picture of quantum mechanics called Tomographic picture [7] , where 
a quantum state is described by a quantity called tomogram. By definition a tomogram 
is a positive measurable probability distribution function of random variables. It has 
been shown that all the physical properties of a quantum system in a given state can 
be studied in terms of its associated tomogram, which then contains the same amount 
of information as the density matrix of the state. 

Although the approach based on tomograms and the one based on density matrices 
are equivalent, it is worth noting two main conceptual differences between these 
two methods: first of all in the tomographic picture one works only with classical 
probability distribution related to relative frequencies of measures. Another difference 
lies in the fact that the construction of a tomogram requires additional information 
about some selected observables in a given experimental setup. 

In the tomographic approach the observables are no longer described by operators: 
they are now associated with functions. The quantization scheme is then based on a 
star product of functions TU]. 

This probability representation of Quantum Mechanics provides us with a useful way 
for analyzing the existence of non local correlations in quantum states of qudits by 
the use of the so called qubit portrait [TT]. This method belongs to the class of 
entanglement-detecting methods based on linear maps. A qubit portrait is indeed 
a linear map defined in the tomographic description of quantum states. Such a 
map constructs a fictions qubit-state associated with the starting qudit-state. This 
construction is needed because Bell's inequality are defined only for systems described 
by dichotomic variables. 

In the following chapters this method will be applied to check Bell's inequalities in 
two classes of quantum states. 



3. Tomographic analysis of entangled states of two qutrits 

We know from Bell's theorem [J] that the existence of nonclassical correlations is given 
by a violation of the inequality 

B<2. (1) 
Using the qubit-portrait method it has been shown [T^] that the Bell's number for the 
maximally entangled state of two qutrits 

I^>-^{|11> + |00> + 1-1-1)} (2) 

can take values greater than 2, for example being B = 1 + y/2. 

In this letter we aim at extending such a result by applying the tomographic 
method to the following two classes of entangled states 

1^^) =^{|l-l) + |00)+e-'^|-ll)}, (3) 

l^^,a) = - 1) + «|00) + e-^'^l 11)}, (4) 

where the real parameters vary in the ranges ip € [0, 27r] and a € [0, 1]. 
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3.1. Ii/'v) state 

In order to apply the tomographic method to a quantum state, first of all we need 
to specify the associated density matrix. For our state, this matrix, expressed in 
the standard basis {|1 - 1), |10), |n), |0 - 1), |00), |01), | - 1 - 1), | - 10), | - 11)}, is 
denoted by and has the following explicit form: 
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(5) 



To obtain the tomogram associated with jl3j . we must evaluate the nine diagonal 
elements of the matrix p^ expressing the density operator p^p after a generic rotation 
of the cartesian reference frame, paramatrized in terms of Euler angles. To this end 
we start by reminding that the matrix rotating the density operator of a single qutrit 
may be cast in the form 



( ie-*("^+x)(l + C0S6') 

^e-*'^ sin 61 
\^ ie'("^-x)(l -cos6i) 



^e-"^sin( 



\/2 

cos Q 
^e^-^sin^ 



ie-»('^-x)(i„cos( 

— ^e*^ sin Q 
ieH'^+x)(i + cos 61 



where 6*, and x denote the three Euler angles [14]. Hence the rotation matrix to be 
used to rotate our two-qutrits state p^p is 

£'p«r(6'l,(/>l,Xi;6'2,02,X2) = D{Bi,(l)x,X\) ® 0(P2Al,X2)- (7) 

where the pedix 1 (2) refers to the qutrit S\ (S2)- Evaluating all the nine diagonal 
elements of the matrix p^ — DpairPv^pair one obtains the following nine-component 
probability vector 

# = {Wi, W2, W3,Wi, W5, We, Wr, Ws, W9) . (8) 

called tomogram of the state p^. We notice that is independent on x due to the 
specific form of p^p. 

To characterize a state through this tomogram amounts at describing it using the 
density-matrix formalism. Nevertheless, if we want to check the Bell's inequality, 
it is easier to start from the tomogram of the state since it allows a systematic 
construction of a new representation mimising a qubit tomogram called qubit portrait 
[11] . In practice starting from our nine-component probability vector we define the 
following new four-component tomogram 

Ijf = {UJI,UJ2,UJ3,^^4), (9) 

where in accordance with 12. we put 

uji^Wi, (10a) 
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UJ2^W2+ W3, (106) 

uj3 = Wi + W7, (10c) 
uji^W5 + We + Ws + Wg. (lOd) 
at represents the tomogram of a quantum state of a fictions two-qubit system. This 
mapping procedure qutrit — >■ qubit provides the way to construct the Bell's number 
of our original state p^p. 

From an experimental point of view to check the Bell's number of a two qubit 
system we must fix for each qubit in a generic way two orientations of our measuring 
apparatus. In this way we perform two measurements for each qubit. We denote by 
a = (^'ia,0ia) and b = (^if,, ^ib) (c = {02c,'t'2c) and d = (92d,4>2d)) the orientations 
of the apparatus corresponding to the measurements performed on qubit 1 (2). To 
construct the expression of B for our state p^p it is useful to introduce the 4x4 matrix 

T = [ UJ^c (11) 
where UJ^ is a 4 x 1 column vector obtained from the vector "at by formally assigning to 



the angle variables 9i, (pi, 62, and (f)2 the generic but fixed values 6 

{^-- 
as 



09 1: and 



a,b;j — c,d). The Bell's number of our original state p^ can then be expressed 



/ = 



(13) 



S-|Tr(/T)|, (12) 
where / is the polarization matrix defined as 

(I -1-1 1 \ 
1-1-11 
1-1-11 

V -1 1 1 -1 / 

It is worthnoting that B is a function of the eight angles 9ia, Ou, 92c, 02d, 4>ia, 4'ib, 
(j)2c-, and (j)2d so that we can look for a violation of Bell's inequality by maximizing 
B with respect to such eight parameters. The result of such a maximization is still 
a function of the state parameter if. Due to the mathematical difficulties of such a 
procedure, we have to resort to a numerical maximization. The function of ip obtained 
in this way is shown in Fig. 1. 



rl 






Figure 1. Bell's number of the state, where ip is on the abscissa the 

numerical maximization, which gives the Bell's number, is on the ordinate. 

It is obvious that the Bell's inequality is violated in the whole range of variability of 
the parameter <p. Hence in our state nonclassical correlations are present no matter 
which value we give to the phase parameter. 
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3.2. \tp^,a) state 

Next, we investigate the existence of nonclassical correlations in the two-parameter 
normahzed state 

\M = -J=^ {|1 - 1) + «|00) + e-'*'| - 11)} , (14) 

where a e [0,1]. We follow the same method used in the previous paragraph 
immediately representing p^^a = |'^i'<^,o)(V'¥>,a| ^ follows: 
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(15) 



Once again we have to rotate this density matrix exploiting Dpair{di, i;^!, Xi; ^2, ?!'2, X2) 

and then taking the diagonal elements of this rotated matrix p^^a = Dl^airPv^a-^pair in 

— y ' ^ ' 

order to construct the tomogram W of our state P(^,o. Proceeding as in the previous 

section, we consider a qubit-portrait of our qutrit state in the same way as we did for 

the \tjj^) state. Using Equations (10), we obtain 



W2 
W3 



W2 

W5 



We 



Ws + Wg. 



Proceeding with the method exposed above, after a numerical maximization on the 
set {Ola, ^16, ^2c) ^2dj (t'la, (f>ib, <p2c, (t>2d} , wc obtain a function of two parameters and 
a. 

Once again an analytical approach is quite complicated and hence after a 
mimcrical treatment we plot this function of both parameters getting the surface 
presented in Fig. 2. 

The most striking feature of this figure is the presence of domains in the {ip, a)-plane 
where a sudden change in the value of Bell's number occurs. In particular, we notice 
the existence of a domain around the point ((po,ao) ~ (Utt/S, 0.9) where the Bell's 
number goes down to a value around 2. In order to see this behaviour in more detail 
we report iu Fig. 3 the Bell's number of our state against a keeping ip fixed at ll7r/6. 
This graph clearly shows the existence of domains in which the Bell's number of 
our state crosses the critical value 2. For example keeping a and <p fixed at values 
(990jQo) = (ll7r/6, 0.9) it is possible to proceed with an exact maximization of the 
Bell's number which gives the result Bq = 1.99999999976005. This means that the 
state satisfies the Bell's inequality in that point. 

This anomalous behaviour may be due to many effects. 

First of all, it could just be an effect of our numerical approach to the problem or 
it could mean that our choice of the qubit portrait is no longer correct in this domain 
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Figure 2. Bell's number for the |V'v>a) state plotted against the two parameters 
ip and a. 
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Figure 3. Bell's number for the {tpip^a) state plotted against the parameter a 
when <p is fixed at the value ipo = ll7r/6. 



of the parameter space, and we should test the value of Bell's number for every other 
possible qubit portraits. However, any other choice in our qubit portrait we tested 
gives a similar result. If this behaviour is not due to numerical maximization problems 
or to unfortunate choices of qubit maps, this could mean that our state undergoes 
a transition from the condition in which nonclassical correlations are present to a 
condition in which the only correlations present in the state are of a classical kind. 
If this is the case then there exists a critical value of the parameter a at which this 
transition takes place. 

It is of interest to observe that the time-dependent state 

\A) - , \ ,^ .{|l-l)+cosM|00)+e-'""/6|-ll)}.(16) 



^2 + cos(w<)2 

is such to span progressively and periodically all the states in the class {|'0ii7r\6,a)i € 
[0, 1]}. As a consequence it evolves from the initial condition 

-zll7r/6| 



100) 



11)}, 



(17) 



in such a way that when ujt ~ arccos(ao) the Bell's number is less than 2. To exploit 
the properties exhibited by \ip, a) thus one might wonder whether an Hamiltonian H{t) 
exists such that \ipit)) satisfies the corresponding Schrodinger equation 

Am) 



dt 



H{t)\m) 



(18) 
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To this end consider the Hamiltonian 

When the coupHng constant \{t) has the form 

^ ^ 5 + cos(2wi) ^ ^ 

and r^o is an arbitrarily assigned energy. 

It is easy to verify that the time evolution operator U{t) associated with H{t) may be 
written down as 

U{t) = e-*""(^^+'^2')*exp (iA(t)[iS'fS'|(S'+S'2-+e-^S'f S'2+)-i(e'^S'+S'2-+S'f S'2+)S'fS'|])(21) 
where K{t) — X{x)dx and that 

u{t)\m) = \m) (22) 

with |^/'(0)) and \'ip{t)) given by (17) and (16) respectively. 

It is not difficult to convince oneself that starting from (19) other examples of Hamil- 
tonians such that \tp{t)) satisfies eq. (18) may be constructed. 



Recently a new method has been proposed to investigate the presence of quantum 
correlations in a given quantum state |15) . This approach suggests to use the non- 
negative quantity Q{A,B) — S{A,B) — I^a^{A,B) to characterize the existence of 
nonclassical correlations in the state of a quantum bipartite system. Here A and B 
are meant to be the two subsystems, >S'(^, B) is the quantum mutual information and 
ImnxiA, B) is the maximum classical mutual information over all possible choices of 
local measurements. The interesting property is that the function Q{A, B) vanishes 
only for states wherein the only correlations are classical. Nevertheless, from a 
practical point of view, the maximization involved in the evaluation of /max(^7 B) 
is often quite complicated. Indeed, classical mutual information /(A, i?), depending 
strongly on the choice of measurement strategies, should be maximized over every 
possible POVM. One route toward the resolution of this problem is to confine 
to projective measures only; this, in turn, allows the evaluation of the quantity 
Q^{A, B) — S{A^ B) — I^^^{A, B) which is meant to surrogate Q{A^ B). Since, indeed, 
the set of projective measures is just a subset of all POVM, we must have 

Q{A,B)<QP{A,B). (23) 

Whether Q^(A, B) may help to extract definitive replies on the existence of quantum 
correlations in a bipartite system in a given quantum state is still an open problem. We 
note indeed that when the function Q^{A, B) is positive, we are not able, in general, to 
give a definitive answer on the presence of quantum correlations since, in this case, we 
do not know anything about the positivity of Q{A, B). Only when Q^{A, B) vanishes, 
we can state with certainty that Q(A,B) = which, in turn, means the existence of 
only classical correlations. 

When applied to the IV'i^.a) state we get Qp{A,B) depending only on a and this 
result is shown in Fig. 4. Since the class \ipip) coincides with \ipip^i), we immediately 
see that Q^{A, B) in this latter case assumes a constant and positive value deducible 
from Fig. 4. 

Hence we see that, for the two-parameter state we are considering, the Q^{A,B) 
function cannot reproduce the result we obtained with our tomographic method when 
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the parameter a crosses its critical value around 0.9. In the case we are considering, the 
method reported in [15j is not able to give any useful answer concerning the presence 
of quantum correlations. 

On the contrary, the tomographic method we applied is a more powerful tool of 
analysis for the state we studied. Indeed, the tomographic approach allows us to give 
definitive replies about the existence of quantum correlations in a wide range of the 
parameter domain where the Bell's number is greater than 2. 

4. Conclusions 

We have examined a class of one-parameter two qutrits state using the tomographic 
approach showing without any doubt the existence of non-classical correlations in all 
the domain of the parameter. 

We have also studied a class of two-parameter states of two qutrits finding almost 
everywhere in the domain of the two parameters non-classical correlations. We 
find in this class also examples of states which might be compatible with the Bell's 
inequalities. This fact suggests the possibility of a passage from states exibiting non- 
classical correlations to states dominated by classical correlations only at the transition 
point ~ (i|^,0.9) when a is varied and ip is kept at its value 

In view of such a property we have demonstrated that the sub-class {|V''iiTr\6.a) i Q G 
[0, 1]} may be dynamically generated and the relevant time dependent Hamiltonian 
accomplishing this target is explicitly given. 

The comparison with another available method [15j for studying nonclassical 
correlations of two-qutrit system shows that the tomographic approach provides a 
tool to detect some peculiarities of the correlations which are not detected by using 
mutual information of bipartite system. The results of this paper can be extended to 
higher spins. 

Acknowledgment s 

V.I.M and M.A.M thank the Russian Foundation for Basic Research for partial support 
under Projects Nos. 07-02-00598 and 09-02-00142 and the University of Palermo (Italy) 
for kind hospitality. 




14 h 



Figure 4. Function QP{A, B) for the I'l/'v'.a) state plotted against the parameter 
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